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SUMMARY 

Concept of using an equivalent monodl sperse spray to represent the vapori- 
zation behavior of polydlsperse sprays has been examined by numerically solving 
two turbulent vaporizing sprays. One Involves the Injection of freon-11 In a 
still environment, whereas the other Is a methanol spray In a still but hot 
environment. The use of three different mean sizes, namely, Sauter mean diam- 
eter, volume median diameter, and surface-area mean diameter, has been Investi- 
gated. Results Indicate a good degree of correlation between the polydlsperse 
spray and Its equivalent monodlsperse sprays represented by the volume median 
diameter and the Sauter mean diameter, the former giving slightly better 
results. The surface-area mean diameter does not provide as good a correlation 
as the other two mean diameters. 


INTRODUCTION 

This work Is aimed at examining the concept of an "equivalent" monodls- 
perse spray with some suitable mean drop size. Purpose Is to explore If an 
"equivalent" monodlsperse spray can adequately represent the behavior of a 
given polydlsperse spray, and If It can, what Is the droplet size of this 
equivalent spray. The work Is Important because the adequate Initial condi- 
tions are currently not available for spray, computations, and because the 
numerical efforts Involved In solving realistic polydlsperse sprays could be 
enormous . 

Most researchers have employed the Sauter mean diameter to represent a 
polydlsperse spray. There Is no evidence, however, to establish that this 
diameter Is the most suitable diameter for simulating the behavior of polydls- 
perse sprays. There have been some studies (refs. 1 to 4) which have attempted 
to find the most representative diameter. Dickinson and Marshall (ref. 1) con- 
cluded that for diffusion-controlled vaporization, no mean diameter can ade- 
quately represent the vaporization characteristics of polydlsperse sprays. 
Alkldas (ref. 2) conducted an analytical study on steady state vaporization 


*Summer Faculty Fellow. Present address: Dept, of Mechanical 

Engineering, University of Illinois, Chicago, Illinois 60680. 



and concluded that the Sauter mean diameter best correlates the overall vapori- 
zation behavior of sprays of different Initial size distributions for both 
diffusion-controlled and radiation-controlled vaporization. A study by 
Aggarwal and Slrlgnano (ref. 3) Indicated that the Ignition characteristics of 
polydlsperse sprays are best represented by the surface-area mean diameter and 
not by the Sauter mean diameter (SMD). A more recent study by Aggarwal 
(ref. 4) concluded that the flame characteristics of polydlsperse sprays are 
best simulated by the SMD, although the surface-area mean diameter Is still 
the more appropriate one for representing the Ignition behavior. Yet another 
Investigation (ref. 5) Indicated that the volume median diameter Is a more 
suitable diameter. Clearly there Is no general consensus on the most optimum 
mean diameter. Moreover, these studies have analyzed Idealized situations. 

The present study is directed toward examining the above Issue by con- 
sidering a more realistic situation of a vaporizing turbulent spray for which 
the Initial droplet size distribution Is obtained experimentally. Two cases 
are considered. First Is an evaporating Freon-11 spray produced by an air- 
atomizing Injector In a still environment. Experimental data for Initial con- 
ditions are available for this case. Second Is an evaporating methanol spray 
which Is also assumed to be produced by an air-atomizing Injector but In a hot 
environment. Since no experimental data Is available for this case, the 
Initial conditions are assumed to be the same as given In the first case. For 
both the cases, results are obtained for the polydlsperse spray and compared 
with those for three equivalent monodisperse sprays represented respectively 
by the Sauter mean diameter (d 3 2 ), the surface-area-mean diameter (d 2 o)t 
and the volume median diameter (do.5v)- 


Brief Description of the Physical Model 

Major assumptions for the continuous phase are: axlsymmetrlc and steady 

flow, boundary layer approximations apply, equal exchange coefficients of all 
species and heat, buoyancy only affects mean flow, and negligible effects of 
mean kinetic energy and radiation. The analysis employs Favre-averaged gov- 
erning equations and a k-c-g turbulent model for the gas-phase, since this 
approach provided good predictions of the structures of variable density 
single-phase jets, as well as evaporating and combusting sprays, during past 
work (refs. 6 and 7). Scalar properties (such as density and temperature) are 
obtained by using conserved scalar approach In conjunction wth state relation- 
ships constructed by adiabatic mixing calculations. 

The liquid phase is treated by solving Lagranglan equations of motion and 
transport for trajectories of a statlstlcally-slgnlficant sample of Individual 
droplets and then computing source terms for mass, momentum and energy due to 
droplets which appear In the governing equations for the gas phase. This 
Involves dividing the droplets Into n groups (defined by position, diameter, 
and velocity) at the Initial condition, and then computing their subsequent 
life-histories In the flow. The stochastlc-separated-f low formulation of 
references 6 and 7 Is adopted for droplet calculations. In this approach, the 
droplets are assumed to Interact with a random distribution of turbulent 
eddies, providing simultaneously for droplet dispersion by turbulence and 
effects of turbulent fluctuations on Interphase transport. The scales and the 
properties of each eddy and the time of Interaction of a droplet with a par- 
ticular eddy are found using the methods described by Shuen et.al. (refs. 6 
and 7). 
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The assumptions for the droplet transport calculations as well as the 
complete formulation of the separated-f low model are given In detail In 
references 6 and 7, and, therefore, will not be repeated here. 


Definition of Mean Diameters 

Various mean diameters for a polydlsperse spray are defined by the 
relation 


X U '{/„ X ' dn /o x3 dn | U <’> 

where dn represents the number of droplets with diameter between x and 
x + dx. 

For the Sauter mean diameter 1 = 3 and j = 2, whereas for the surface- 
area mean diameter 1 = 2 and j = 0 . 

Normalized cumulative volume for the liquid phase Is defined as 


C v (x) = ) x 3 dn/f x 3 dn (2) 

J Q • J 0 

and the volume median diameter Is defined as the value of upper limit In the 
above Integral corresponding to C v = 0.5. 


Initial Conditions for the Monodlsperse Sprays 

Using the experimental size distribution, the values of d 32 and d 2 o 
were calculated as a function of jet radial location. The result Indicated 
that d 32 and d 20 are almost uniform In the radial direction Initially. 

Thus a single droplet size Is sufficient for specifying the Initial liquid mass 
distribution. 

The volume median diameter can be calculated by using equation ( 2 ), where 
the Integral Is computed numerically by using the experimental size distribu- 
tion. The calculated values for d 2 Q. CI 32 * ar) d ^0.5v f° r the experimental 
size distribution are 48.4, 59.8, and 65.0 pm, respectively. For the monodls- 
perse sprays, the droplet velocity distribution Is obtained by using the 
assumption that the velocity distribution for a given drop size Is the same as 
the experimental velocity distribution for that size. For example, for mono 
disperse spray with d 32 (= 59.8 ym), the Initial velocity distribution 
corresponds to the experiment values for droplet size of 59.8 ym. A slightly 
better velocity distribution may be obtained from the total liquid-momentum 
flux considerations, I.e., by equating the total momentum flux of the monodls 
perse spray to that of the polydlsperse spray. However, the differences 
between the two distributions were found to be small. 
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Results for the Freon-11 Spray 


For the evaporating freon-11 spray, the Initial conditions for two-phase 
calculations are taken from the experimental data of reference 6 at x/D = 50. 
Results for the polydlsperse and three equivalent monodisperse sprays repre- 
sented by d32, d20» and do.5v are presented In figures 1 to 3. Variation 
of freon vapor mass fraction In the jet Is given In figure 1. Agreement 
between the four cases Is excellent, both In terms of the axial and radial 
variations. The profiles of gas temperature In the axial and radial directions 
are also almost Identical for the four cases (not shown). Such an excellent 
agreement may be somewhat misleading because the freon vapor mass fraction at 
the Initial station (x/D = 50) Is relatively high. Its variation Is then con- 
trolled by the fluid-dynamical and mixing processes rather by the liquid-phase 
processes. Figures 2 and 3 are perhaps more representative of the differences 
between the polydlsperse and the equivalent monodisperse sprays. Radial pro- 
files of total freon mass fraction at x/D = 400 and 500 are shown In figure 2. 
Differences, which are relatively large at x/D = 500 are somewhat exaggerated 
because 80 percent of liquid mass Is evaporated and only very large droplets 
are remaining at this location. Thus, figure 2(b), which Indicates a good 
degree of correlation between the polydlsperse spray and Its equivalent mono- 
disperse sprays, provides a better representation of the differences In most 
of the jet. Radial variation of liquid flux at two axial locations Is given 
In figure 3. Again, figure 3(b), which provides a better representation of the 
differences than figure 3(a), shows that the vaporization behavior of polydls- 
perse sprays can be reasonably represented by the equivalent monodisperse 
sprays. Amongst the three monodisperse sprays, the one with the volume median 
diameter seems to provide the closest predictions to the polydlsperse spray, 
whereas the d 2 g sprays shows the largest differences. Variation of liquid 
flux along the jet axis, given In table I, demonstrates similar degree of cor- 
relation between polydlsperse and three monodisperse sprays, with do.5v and 
d20 sprays showing, respectively, the best and the worst agreement with the 
polydlsperse spray. At x/D = 70, d0.5v and d32 sprays overpredict the 
centerline liquid flux. This Is perhaps the fluid-dynamic effect, as the jet 
Interior Is relatively cold and the vaporization rate Is small. Smaller drop- 
lets tend to move away from the centerline whereby reducing the liquid mass 
flux there for the polydlsperse case. For x/D between 100 and 400, the 
do.5v spray continues to overpredict the centerline liquid flux, which Is now 
due to the vaporization effect, l.e., the smaller droplets In the polydlsperse 
spray are vaporizing faster. For x/D greater than 500, the liquid flux Is 
higher for the polydlsperse case because only the large droplets are remaining. 
Flowever, the differences are not Important there, as most of the liquid mass 
has already vaporized. 

Differences between the polydlsperse and d 2 Q sprays are small 
Initially. Further downstream, x/D between 150 and 250, the d20 spray over- 
predicts the centerline liquid flux due to the faster vaporization of smaller 
droplets In the polydlsperse case. For x/D greater than 250, liquid flux Is 
higher for the polydlsperse spray for the same reason as mentioned earlier. 
Perhaps, the radial distribution of liquid flux, as given In table 1(b) at 
x/D = 150, provides a better representation of the differences between the 
four cases. In the Inner region, the monodisperse sprays generally over- 
predict, whereas In the outer region they underpredict the liquid flux as com- 
pared to the polydlsperse case. The differences are due to the vaporization 
effect and the radial distribution of gas temperature. The Inner region is 
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relatively cold, thus, only the smaller droplets In the polydlsperse spray 
experience significant vaporization. In the outer region, the vaporization 
rate Is higher for the monodlsperse sprays as compared to the polydlsperse 
spray since the gas temperature Is relatively higher. 

Total freon mass fraction values at the centerline for the four cases are 
shown In table 1(a). Obviously the differences between polydlsperse and mono- 
dlsperse sprays are much smaller as compared to those In the liquid flux 
comparison. 


CONCLUDING REMARK 

The major conclusion from the above results Is that the vaporization 
behavior of polydlsperse sprays can be represented by the equivalent monodls- 
perse sprays. There Is a reasonably good correlation, both In terms of gas- 
phase as well as liquid-phase properties, between the polydlsperse and three 
monodlsperse sprays. Moreover, the use of equivalent monodlsperse spray 
reduces the computational effort by a factor of six. For the present results, 
the use of volume median diameter best represents the polydlsperse spray 
behavior, although the Sauter mean diameter and the surface-area-mean diameter 
also provide a reasonably good correlation to the polydlsperse spray. It 
should be noted, however, that the freon-11 spray, due to Its fast vaporization 
characteristics and due to high value of Initial fuel vapor concentration, may 
not be a good test case for examining the degree of correlation between the 
polydlsperse spray and Its equivalent monodlsperse sprays. This Issue Is pur- 
sued further by using a methanol spray, as discussed In the next section. 


Methanol Spray Results 


Results are now presented for a vaporizing turbulent methanol spray. The 
purpose Is to further examine the degree of correlation between the polydls- 
perse spray and Its equivalent monodlsperse sprays. The physical model Is 
essentially the same as that for the freon-11 case, except that a methanol 
spray In hot still surroundings, which Is a temperature of 800 °K, Is consid- 
ered. Initial distribution of gas-phase properties Is assumed to be the same 
as that for the freon-11 spray. Initial droplet size and velocity distribu- 


tions are also assumed 
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no other experimental data Is available in the literature. Another reason for 
using the same Initial conditions as much as possible Is to alter only the 
vaporization characteristics of the spray without changing Its fluid- dynamical 
characteristics. Methodology for calculating droplet size-history and trajec- 
tory Is Identical to the one employed In reference 7. Transport and thermo 
dynamic properties as well as the equations of states (for calculating gas 
temperature and density) are also computed by following the procedure of 
reference 7. Direct droplet Interaction Is neglected and the thin skin approx- 
imation (ref. 7) Is employed. Further details are given In the above 
reference. 


Results are again obtained for the polydlsperse case and compared with 
those for the three monodlsperse sprays. Methanol vapor distribution In the 
jet for all four cases Is shown In figure 4. The differences In the radial 
distribution of methanol vapor are negligible. In the axial direction, the 
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agreement between the polydisperse, d 32 , and dg 5 V sprays Is quite reason- 
able. However, the d20 spray Indicates significant departure from the poly- 
disperse spray, the maximum difference being 35 percent as compared to 
10 percent for the d32 and dg. 5 V sprays. It is also interesting to note 
that the differences between polydisperse and monodlsperse sprays are now more 
significant as compared to the freon case. Compare figures 1(a) and 4(a). 

This Is due to the large amount of Initial fuel vapor In the freon case; the 
methanol vapor mass fraction at the Initial location Is zero. Consequently, 
the differences In the vaporization behavior of polydisperse and monodlsperse 
sprays are better highlighted with the methanol spray. More Interesting obser- 
vation Is that even for the methanol case, the results for d32 and dg. 5 V 
sprays are quite encouraging. 

Radial variation of total methanol mass fraction and of liquid flux are 
shown In figures 5 and 6 , respectively. 

Radial distribution of liquid flux at two axial locations Is given In 
figure 5. At x/D = 150, the dg 2 predictions seem to be the closest to 
the polydisperse results. Except near the jet edge, the liquid flux values are 
higher for the dg 2 and dg 5 V sprays because the small droplets In the 
polydisperse spray are vaporizing faster. Note that, unlike the freon case, 
the Interior of the jet Is hot and the vaporization there Is significant In the 
present case. Moreover, the effect of gas-phase convection on droplet vapor- 
ization is the largest at the jet axis. As a result, the liquid flux at the 
axis is lower as compared to the values In the jet Interior. The liquid flux 
values for the d 2 g spray are underpredicted, which Indicates faster vapor- 
ization for this case as compared to the polydisperse case. It is also note- 
worthy that the actual differences between the polydisperse and the monodls- 
perse sprays are smaller than those Indicated In figure 5(a). This Is because 
the liquid flux has been normalized by the center-line value, which Is smaller 
for the dg 5 V and d 32 sprays but higher for the d 2 g spray as 
compared to the polydisperse case. The center-line liquid flux is given In 
table II. 

The radial distributions of methanol liquid flux at x/D = 300 (fig. 5(b)) 
Is remarkedly different from those at x/D = 150. Results for all three mono- 
dlsperse sprays are now in agreement with each other, but are consistently 
underpredicted as compared to the polydisperse case. The underprediction is 
due to the fact that most of the liquid (more than 80 percent) has already 
vaporized at this location and only a few large droplets In the polydisperse 
spray are remaining. Thus, the results In figure 5(a) along with table II are 
more representative of the differences between the four sprays as compared to 
those In figure 5(b) . 

Figure 6 shows the radial variation of total methanol mass fraction. 
Obviously, the agreement Is much better here as compared to that In figure 5. 
Table III, which gives the axial variation of total liquid and vapor mass flow 
rate, provides a global comparison of the vaporization behavior of polydisperse 
spray with Its equivalent monodlsperse sprays. Again, the degree of correla- 
tion between the polydisperse and dg 5 V sprays, and between the polydis- 
perse and d32 sprays Is quite acceptable up to x/D = 250. It deteriorates 
further downstream but is unimportant as much of the liquid has already vapor- 
ized. Another interesting observation emerges from the comparison of figures 3 
and 5. The differences between the polydisperse and monodlsperse sprays are 
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relatively larger for the methanol case and are more typical of the degree of 
correlation between the vaporization behavior of polydlsperse spray and Its 
equivalent monodlsperse spray. The methanol results are nevertheless encour- 
aging, specially for the dg. 5 v and CI 32 sprays. 


CONCLUSIONS 

Structure of evaporating turbulent sprays Is numerically computed and the 
use of “equivalent" monodlsperse spray for simulating the vaporization behavior 
of realistic polydlsperse sprays Is Investigated. The physical and numerical 
models are the same as In reference 6 . Three mean diameters which are examined 
are the Sauter mean diameter, volume median diameter, and surface- area mean 
diameters. Results are obtained for a freon-11 spray In a still environment 
and for a methanol spray also In a still but hot environment. For both cases, 
the results for the polydlsperse spray are compared with those for Its equiva- 
lent monodlsperse sprays. Major conclusions are: 

(1) For the freon-11 case, the degree of correlation between the polydls- 
perse and equivalent monodlsperse sprays Is quite good. Amongst the three mean 
diameters, the volume median diameter and the surface-area mean diameter pro- 
vide respectively the best and the worst representation of polydlsperse spray 
behavior. The results for the Sauter mean diameter are quite close to those 
for the volume median diameter. 

(2) For the methanol case, the correlation between the equivalent mono- 
dlsperse and polydlsperse sprays Is not quite as good as for the previous case 
but Is still acceptable. Again, the volume median diameter best simulates the 
vaporization behavior of polydlsperse spray. 

(3) The methanol spray results are perhaps more typical of a turbulent 
vaporizing spray than those for the freon 11 because of the high Initial vapor 
concentration In the latter. 

(4) The use of a suitable mean diameter Is also encouraging due to the 
reduction In computational efforts by a factor of six. This will become more 
of an Issue as advanced vaporization models (ref. 8 ) are employed or If the 
combustion situation Is considered. 

(5) The general conclusion Is that the use of a suitable defined monodls- 
perse spray for representing a realistic polydlsperse spray Is quite encour- 
aging. The degree of correlation In the more Interesting region of spray, 
where more than 80 percent of the liquid vaporizes. Is acceptable to quite 
good. It can perhaps be Improved further by employing two mean diameters. 
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TABLE I 


(a) Variation of liquid flux and total freon mass fraction 


along the jet axis 


x/D 

Liquid flux at 
the jet axis 

(normalized by the value at 

x/D = 50) 

Total freon mass 
fraction at the jet axis 


Poly- 

disperse 

d 0.5v 

D 

a 

Poly- 

disperse 

m 

m 

d 20 

m 

1.0 

l 

1 

1 .0 


0.526 

0.526 

0.526 

EH 

.560 

umSm 

1 

.560 

.385 

.390 

.390 

.386 

KB9 

.328 

■ 

.318 

.314 

.287 

.285 

.285 

.289 

150 

.194 

Hfl 

.160 

.287 

.216 

.214 

.204 

.250 


.102 

.m 


.141 

.165 

.167 

.191 

.203 


.032 

.034 


.019 

.119 

.124 

.126 

.124 

500 

.012 

wii 


.001 

.095 


.096 

.09 

600 

.005 

■H 


0 


.078 

.075 

.073 


(b) Radial variation of liquid 
flux at x/D = 150 


r/x 

Polydlsperse 

d 20 

d 0.5v 

d 32 

0.013 

1 



0.160 

.044 

i ■ 

.218 

.227 

.27 

.082 

.10 

.12 

.12 

.09 

.12 

.03 


.034 

.02 

.17 

.003 

mm 

.001 

.001 


TABLE II. - LIQUID FLUX AT THE CENTERLINE, 
NORMALIZED BY THE VALUE AT x/D = 50 


x/D 

Polydl sperse 

d 0.5v 

d 32 

d 20 

50 

1 .0 

1 .0 

1 .0 

1.0 

70 

. 539 

cco 

551 

coo 

100 

.295 

.275 

.245 

.273 

150 

.159 

.133 1 

.137 

.225 

200 

.101 

.12 1 

.116 

.119 

250 

.056 

.086 

.078 

.045 

300 ! 

.029 

.038 ' 

.026 

.006 

350 

.015 

.012 

.005 

.0003 
































# 


TABLE III. - COMPARISON OF TOTAL LIQUID 
MASS FLOW RATE AND TOTAL VAPOR MASS 
FLOW RATE FOR THE FOUR SPRAYS 


(a) Total liquid mass flow rate 


x/D 

Polydl s- 
perse 

d 32 

d 0.5v 

d 20 

Total liquid mass flow rate, kg/sec 

50 

2.39x10-4 

2.39x10-4 

2.39x10-4 

2.39x10-4 

70 

2.09 

2.04 

2.12 

1.89 

100 

1.77 

1.74 

1 .83 

1.5 

150 

1 .34 

1.27 

1.38 

.98 

200 

.98 

.86 

1 .0 

.58 

250 

.70 

.52 

.66 

.27 

300 

.47 

.26 

.38 

.07 

350 

.30 

.09 

.17 



(b) Total 

vapor mass flow rate 

x/D 




d 20 

Total vapor 

mass flow 

rate, kg/sec 

50 

0 

0 

0 

0 

70 

.25x10-4 

.22x10-4 

.21x10-4 

.31x10-4 

100 

.56 

.51 

.49 

.69 

150 

.99 

.99 

0.94 

1 .21 

200 

1.34 

1 .39 

1.33 

1.61 

250 

1 .62 

1.73 

1.67 

1 .92 

300 

1.85 

1 .99 

1 .95 

2.1 

350 

2.01 

2.15 

2.15 

2.15 
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(b) X/D = TO. 

Figure 3.- Radial variation of mean liquid freon 
flux for the polydisperse and equivalent mono- 
disperse SPRAYS. 
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(b) Radial variation at x/D = 250. 

Figure 4.- Variation of methanol vapor mass frac- 
tion in the turbulent evaporating spray for the 
polydisperse and three equivalent monodisperse 
sprays. 


LIQUID FLUX, G/G. 




(B) x/D = 300. 

Figure 5.- Radial variation of mean liquid metha- 
nol flux for the polydisperse and equivalent mono- 
disperse SPRAYS. 



Figure 6.- Radial variation of total methanol mass 

FRACTION FOR THE POLYDISPERSE AND EQUIVALENT 
MONODI SPERSE SPRAYS, X/D = 300. 
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